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Abstract 

We investigate the electromagnetic field generated by a point charge moving along a 
helical trajectory inside a circular waveguide with conducting walls filled by homogeneous 
dielectric. The parts corresponding to the radiation field are separated and the formulae 
for the radiation intensity are derived for both TE and TM waves. It is shown that the 
main part of the radiated quanta is emitted in the form of the TE waves. Various limiting 
cases are considered. The results of the numerical calculations show that the insertion of the 
waveguide provides an additional mechanism for tuning the characteristics of the emitted 
radiation by choosing the parameters of the waveguide and filling medium. 



PACS number(s): 41.60.Ap, 41.60.Bq 



1 Introduction 

A charged particle confined to the helical orbit is a source of high intensity electromagnetic 
radiation over a broad range of wavelengths with a number of remarkable properties such as high 
collimation and high degree of polarization (see, for instance Refs. [UG1E]). These properties 
have resulted in extensive applications of this radiation in a wide variety of experiments and in 
many disciplines. In particular, the helical motion of an electron beam is employed in helical 
undulators to produce circularly polarized radiation in narrow angular cone in the forward 
direction (5J El El [8] . In the proposal of Ref . [4] the helical undulator radiation was used to 
generate a polarized positron beam. Synchrotron radiation from relativistic electrons spiralling 
in magnetic fields is the main mechanism to explain the emissions of many objects in radio 
astronomy (see [9] and references therein). Most of the works on the radiation from the helical 
trajectory refer to radiation in free space. It is well known that the presence of medium can 
essentially change the characteristics of the electromagnetic processes and gives rise to new types 
of phenomena such as the Cherenkov, transition, and diffraction radiations. In particular, the 
operation of a number of devices assigned to production of electromagnetic radiation is based 
on the interaction of charged particles with materials (see, for example, |10j). 

The synchrotron radiation from a charged particle circulating in a homogeneous medium 
was considered in Ref. [llj. In this paper it was shown that the interference between the 
synchrotron and Cherenkov radiations leads to interesting effects. New interesting features arise 
in inhomogeneous media. In particular, the interfaces of media can be used to control the 
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radiation flow emitted by various systems. In a series of papers started in Refs. |12[ |13j . we 
have considered the simplest geometries of boundaries between two dielectrics with spherical and 
cylindrical symmetries. The synchrotron radiation from a charge rotating around a dielectric 
ball enclosed by a homogeneous medium is investigated in Refs. |13|. I14j . It was shown that 
if for the material of the ball and the particle velocity the Cherenkov condition is satisfied, 
strong narrow peaks appear in the radiation intensity. At these peaks the radiated energy 
exceeds the corresponding quantity in a homogeneous medium by several orders of magnitude. 
A similar problem with the cylindrical symmetry has been discussed in Refs. |12[ \15\ [161 117j . 
In Ref. [12] we have developed a recurrent scheme for constructing the Green function of 
the electromagnetic field for a medium consisting of an arbitrary number of coaxial cylindrical 
layers. The investigation of the radiation from a charged particle circulating around a dielectric 
cylinder immersed in a homogeneous medium, has shown that under the Cherenkov condition 
for the material of the cylinder and the velocity of the particle, there are narrow peaks in the 
angular distribution of the number of quanta emitted into the exterior space. For some values 
of the parameters the density of the number of quanta in these peaks exceeds the corresponding 
quantity for the radiation in vacuum by several orders. The radiation by a longitudinal charged 
oscillator moving with a constant drift velocity along the axis of a dielectric cylinder immersed 
in a homogeneous medium is investigated in Refs. [18, 19J. As in the case of the circular motion 
it was shown that the presence of the cylinder provides a possibility for an essential enhancement 
of the radiation intensity. The properties of the radiation from a charged particle moving along 
a helical orbit in homogeneous dispersive medium are investigated in [20]. The corresponding 
problem for the charge moving in vacuum has been widely discussed in literature (see, e.g., Refs. 
[H El El [21] and references given therein) . The electromagnetic field and the radiation in the case 
of particle following the helical path inside a dielectric cylinder immersed into a homogeneous 
medium are studied in Refs. |22[ 123] . Recently the influence of a homogeneous transparent 
medium on the radiation of relativistic particles in planar undulators is considered in Ref. [24] . 

In the present paper, we study the electromagnetic field and the radiation intensity for a 
charge moving in a helical orbit inside a circular waveguide with dielectric filling. Note that 
the radiation parts of the fields in the case of vacuum inside the waveguide are investigated 
in Ref. [25]. The plan of the paper is as follows. In section [2] we derive expressions for the 
electric and magnetic fields by making use of the corresponding formulae from [23] for the 
geometry of a dielectric cylinder immersed into a homogeneous medium. Analytic properties 
of the corresponding Fourier components are investigated. In section [3] the radiation fields 
are separated and they are presented as a superposition of the waveguide eigenmodes. The 
formulae are derived for the radiation intensity of TE and TM waves and numerical examples 
are presented. Section S] concludes the main results of the paper. 



2 Electromagnetic fields inside a waveguide 

Consider a point charge q moving along the helical trajectory of radius po inside a circular waveg- 
uide with conducting walls. We will denote by p% the radius of the waveguide and will assume 
that it is filled by homogeneous dielectric with permittivity Eq. The particle velocities along 
the axis of the waveguide (drift velocity) and in the perpendicular plane we denote by vn and 
v±, respectively. In a properly chosen cylindrical coordinate system (p,cp,z) the corresponding 
motion is described by the coordinates 

P = P0, (f) = LU t, z = v\\t, (1) 
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where the z-axis coincides with the waveguide axis and loq = v±/pq is the angular velocity of 
the charge. This type of motion can be produced by a uniform constant magnetic field directed 
along the axis of a cylinder, by a circularly polarized plane wave, or by a spatially periodic 
transverse magnetic field of constant absolute value and a direction that rotates as a function 
of the coordinate z. In the helical undulators the last geometry is used. 

In accordance with the symmetry of the problem, we present the electric and magnetic fields 
in the form of the Fourier expansion 

F,(r,t) = Y, j m{ *-" ot) / dk z e tk ^~ v ^F ml (k z ,p) 
m=-oo J - 00 

= 2Re^ e M4>-"ot) / dk z e ik ^ z - v \^F ml {k z ,p), (2) 

where I = p,(p, z, and in the discussion below F = E,H for the electric and magnetic fields, 
respectively. The expressions for the Fourier transforms F m i(k z ,p) are obtained from the corre- 
sponding formulae derived in Ref. [23] for the geometry of a dielectric cylinder with permittivity 
£q immersed into a medium with permittivity E\, taking the limit e\ —* 00. As a result of this 
limiting procedure, the Fourier transforms are written in the decomposed form 

F ml (k z ,p) = F^(k z ,p) + F^(k z ,p), (3) 

where the part F^(k z ,p) corresponds to the fields generated by the charge in a homogeneous 
medium with permittivity Eq and the part F^j (k z , p) is induced by the presence of the waveguide. 
In the case p < po for the homogeneous part of the magnetic field one has 

< ] i = £ P-- 1 ^ WAop), 1 = PA, (4) 



2m a l 

<?A C 



p=±i 



ffffl = E PD { yj m (Xo P ), (5) 

p=±i 



with the coefficients 



n(Qp) = I_ 
2ic 



v±H m+p (X p ) - vn^-H m (X po) 

K z 



(6) 



In these expressions a p = 1, = 2, J m (x) is the Bessel function, H m (x) = H^\x) is the 
Hankel function of the first kind, and 

uj 2 (k z ) 

Aq = m 2 2 £0 - k 2 z , u m (k z ) = mw + k z v\\. (7) 



The part H^j(k z , p) induced by the presence of the waveguide is given by the formulae 



The corresponding expressions for p < pq are obtained from ([4]), ([5]) by the replacements J T^- H. 
r (i) _ qk z 

p=±l 



H% = E^" 1jD -WAop), 1 = PA, (8) 

H ™ = E pD^Jm(Xop), (9) 

p=±i 
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where we have introduced the notation 

Aq H m (Xopi) 



= — 



7T 

2ic 



k z Jm(AoPl) 
iv±p J m -p(A pi) 



Jrni^OPo) — V±J m +p{^OPo)' 



H m + p (XoPl) 



Ti 



vrpiAo Jm(AoPi) J^(Aopi) J m+ z(A /9i 



(10) 



and the prime stands for the derivative with respect to the argument of the function. 

The part in the electric field due to the presence of the waveguide is given by the formulae 



E 



(i) 

ml 



qci 



l-eri 



4iru m (k z )eo 

qicXpk z 
2iroj m (k z )e 



^2 p ai J m+p (X p) 



p=±i 



u 2 m {k z )e Q 



+ kl)D$-XlD, 



(-p) 

m 



D$J m (Xop), 



11) 



(12) 



p=±i 



where I = p,<f>. The corresponding formulae for the part E^j are obtained from pip . (|12p by 

the replacement Dm — > Dm ■ Note that in the limit po — ► Pi one has — ► —Dm P ^ and the 
fields vanish. We could expect this result, as when the charge is on the surface of the waveguide 
the charge and its image compensate each other. 

The formulae given above describe the total electromagnetic field inside the waveguide. To 
separate the parts corresponding to the radiation we need the analytic properties of the Fourier 
components as functions on k z . From formula (|10p it follows that the function D$ has singular 
points corresponding to the zeros of the function J m (Xopi) and its derivative. Note that in the 
second and third summands on the right of formula (|10p the singularities at the zeros of the 
functions J m ±i(Ao/Oi) cancel out and the function D m is analytic at these points. We denote by 
jmjn, n = 1, 2, . . ., the nth positive zero of the Bessel function (a = 0) and its derivative (a = 1): 



j£)(X oPl ) = J^ti$ n ) = 0, 



(13) 



where Jm{%) = d a J m (x) / dx a . These zeros describe the eigenmodes for the cylindrical waveg- 
uide and are known as TM modes for the case a = and as TE modes for the case a = 1 
[26]. The corresponding modes for the projection of the wave vector on the cylinder axis are 
determined from the relation AqPi = jm,n by taking into account expressions ([7]): 



£o(3\\k z + m 



2 W £ 



This equation has real solutions under the condition 

&(*)2 (i-rA < 1 />M : 
u m,n I 1 h>\\ ) L i u m,n ■ 



Jm,n 



CJm,n 



0, 



£0- 



(14) 



and these solutions have the form 



u — i.O,±) 



muW £ o 



c 1 



(15) 



(16) 



The values km',n ^ correspond to the simple poles of the functions E m i(k z , p) and H m i{k z , p). Note 
that if the Cherenkov condition for the velocity of the particle along the axis of the waveguide 
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is satisfied, j3\\ > 1, then inequality (|15p is valid for all values n = 1,2, In the case (3\\ < 1, 

condition (|15p determines the maximal value for n, which we will denote by ni 



'max- 



. ((7 ) m^ppi^/go .(a) 
^ M < / t < 3 ( CT ) • 

max T J- 



The value of the function uj m (k z ) at the points fcm'n^ is equal to 



(17) 



(18) 



For m = 0, equation (|14p has real solutions only when /3m > 1 and 



k z — k, 



0,71 



> W 0,n 



Wife, 



(«i±) 
0,n • 



(19) 



For real values Eq an d under the condition (|15p the poles are situated on the real axis of the 
complex plane k z . In formula ([2j) it is necessary to give the way by which these poles should 
be circled. For this we note that in physical situations the dielectric permittivity is a complex 
quantity, eq = e' + ie^, and the imaginary part e'q determines the absorbtion in the medium. 
Under the condition |e '| <C £q, for the imaginary part of km',n from (|16|) one has 



Imfc&£> =±Cie[, , (u; m ), C x > 0, 



(20) 



where e'q (u ro ) ^ for u> m ^ 0. It can be also seen that (bm,n — l)Re fcm,'n > 0. If the Cherenkov 
condition is not satisfied, = Unyi^/c < 1, from the formulae given above it follows that 



Re > 0, u; m (Re />;±) ) > 0, ±Im k^ ] > 0. (21) 

In this case the poles (km\r?) are situated in the upper (lower) half of the complex plane 

k z . In the limit Eq — > 0, deforming the integration contour we obtain the rule for avoiding the 
poles plotted in figure [TJ In the case /3m > 1 the corresponding inequalities have the form 

ReAfc^ < 0, Tw m (Refc(gO) < o, ImA;^ < 0. (22) 

In the way similar to that for the previous case, deforming the contour for the integration over 
k z , we obtain the avoidance rule for the poles given in figure [2j 



3 Radiation intensity in the waveguide 
3.1 Radiation fields 

In this section we consider the radiation field propagating inside a cylindrical waveguide with 
perfectly conducting walls for large distances from the charge. First of all let us show that 
the part corresponding to the first term on the right of formula ([3]) does not contribute to the 
radiation field. This directly follows from the estimate of the integral over k z on the base of the 
stationary phase method. As in the integral over k z the phase k z z has no stationary points, for 
large values \z\ the integral vanishes faster than any degree of l/|z|, under the condition that the 
preexponential function belongs to the class C°°(R). It follows from here that the radiation field 
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Figure 1: Integration contour in the k z plane for the case /3n < 1. 
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Figure 2: Integration contour in the A; 2 plane for the case /3n > 1. 



is determined by the singularity points of the preexponential function. As it has been mentioned 

(p) 

before, for the integral over k z the only singular points are the poles of the function D m at the 

>.±) 



points k z = km,n , determined by relations (|T6|) . To find the corresponding contribution, we 
note that the integration contour over k z has the form depicted in figure [T] for (3\\ < 1 and the 
form depicted in figure [2] for the case /3\\ > 1. At large distances from the charge the integration 
contour can be closed by large semicircle in the upper (lower) half-plane for z > v\\t (z < i>ii£). 
This choice is caused by the fact that the integrand exponentially vanishes in the upper (lower) 
half-plane. As a result when the radius of the large semicircle goes to infinity the corresponding 
integral vanishes. Hence, for large values \z\ the integral over k z , by residue theorem, is equal to 
the sum of residues inside the contour multiplied by 27risgn (z — vut). For e'q — ► and large n, 

when b m ] n > - 1 , the poles kt,t ] have a finite imaginary part and the corresponding 

contribution exponentially vanishes in the limit z — > oo. As a result these poles do not contribute 
to the radiation field. 

As before we consider two cases. When /3n < 1 we close the integration contour in figured] 
by the semicircle of large radius in the upper half-plane for z > vnt and in the lower half-plane 
for z < vnt. As a result for the radiation field one finds 



where a = + (a 



Fi(r,t) = aRe 



00 "mix 



<x=0,l m=0 n=l 

) corresponds to the case z > v»t (z < vnt) and we use the notation 

F m T ; a) (r,t)=4ni Res F ml (k z , p) e *(w+*.*-«mt) . 
i. _i.( CT > a ) 



(23) 



(24) 
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Each term in the sum on the right of formula (|23p describes waves with the frequency u>m',n 
propagating along the positive direction of the axis z for a = + and for a = — , 1 < &m,« < 
l/^/l — and waves propagating along the negative direction of the axis z for a = — , &m,n < 1- 
If the Cerenkov condition is satisfied, /3u > 1, then closing the integration contour in figure [2] by 
large semicircle in the upper or lower half-plane in dependence of the sign for z — v»t, for the 
vector potential of the radiation field one finds 



Fi(r,t) = -6{v\\t - z)Re 



00 "mix 



E E EttfM) 

a=±<T=0,lm=0 n=l 



(25) 



where 9{x) is the Heaviside unit step function. Separate terms in the sum in formula (|25p 
describe waves with the frequency l^^n | propagating along the positive direction of the axis 
z for a = + and for a = — , bm]n > 1, and along the negative direction of the axis z for a = — , 
b { m\ < 1. Note that for 

frm,n > 1 we have no waves propagating along the negative direction of 
the axis z. By taking into account the formulae for the Fourier components of the fields and 
evaluating the residues by the standard formulae of the complex analysis, we find the following 
expressions for the radiation parts of the z-components of the fields: 

Et a) (r,t) = 4gJ ™ ( f nPQ lf Jmi&p/Pi) ex P [^(m y + t®$z - agt)j, (26) 

H mz [r,t) - « (1 ) (1) Jm{3 m , n p/Pl) 

\ /t 0Pl\Jm,n III- )-J m {Jm,n) 



exp[i(rrup + km$ z - uj&ftt)] 

I ) K^'j 

mu Jl + b$n 

and Hmf ^ (r, t) = E^z^ (r, t) = 0. The transverse components are found from the formulae 

E^(r,t) = ^^ (0 \H^(r,t) = £^e 3 xE^(r,t) ; (28) 

Jm,n CKtf 



H^f } (r,t) = ^^«,ESf ) (r,t) = -^||e 3 xHSf ) (r,t) J (29) 



i iii.n 



where V (0) = E^f\r,t) and = H&f\r,t), V t = (d/dp,im/p,0), and e 3 is the unit vector 
along the axis of the waveguide. 

3.2 Radiation intensity 

As we have seen, inside the waveguide the radiation is presented in the form of waves with 
discrete set of the values for the projection of the wave vector on the waveguide axis, kz — km,n 5 
n = 1,2,..., determined by formula (|16p . Having the radiation fields we consider the mean 
energy lost per unit time, 



^ I ti>P E 'P+jz E z)P d pd<pdz, T = 2tt/uj . 



(30) 
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The radiation intensity / is presented as the sum of intensities on separate modes: 

oo nl CT) 



I ~ XT ^™ ' ^™ ~ XT ^m',n • (31) 



m=0 <x=0,l a=± n=l 

The term with m = is present only then the condition /Jy > 1 is satisfied and the corresponding 
parts have the form 

, 0) _ 2g\ JHCp,,/Pi) t(1 , _ 2qV 1 v ll JjigjpM 
Each term in the sums of these formulae corresponds to the radiation with the frequency given 



by the formula v\\jon/(piy//3^ — 1). For the radiation intensities on harmonics m 7^ one has 
the formulae 

o 2 t2 ( -(0) / \ i (0,a)| / 

r(0,o) _ 2 ? c Jm(Jm,nPo/Pl) \U m , n \ ,,.(0)2^ 

Vn " 3/2 2 ,(0)2 2 , (0) , mW() V i + W ^ W 

9 „2„,2 -(1)2 .'2/ •(!) / n \ \i 
j(l,a) Z( l V ±_ Jm,nJ m {]m,nP$/ Pi) \^m,n 



c ^pI (iS - m^) j* ($„) ^i + ^g^ 



i 



(34) 



For /3|| < 1 the upper limit of the summation over n is defined by relation (|17p . Otherwise this 
limit is determined by the dispersion law for the dielectric permittivity eq through the condition 
v\\y/eo > c. It can be seen that for the case Sq = 1 the expression for ^ a Imffl obtained from 
formulae ([33]) and (f34"j) coincides with the corresponding formulae in Ref. [25]. Taking /3\\ = 0, 
from formulae (|33|) and (|34p we obtain the corresponding results for the radiation from a particle 
circulating in the plane perpendicular to the waveguide axis {281 . 

In accordance with (|15p . for given m and n the necessary condition for the presence of the 
radiation is the condition 

rrnooPi V^/c > j$ n . /l - /3 2 . (35) 



Now by taking into account the relation jinn ^ tj. for the zeros of the Bessel functions, we 



conclude that under the conditions /?y < 1 and j3±_ < — /3 2 po//?l> with = v±^/eq/c, there 
is no radiation inside the waveguide though the particle moves with acceleration. If /3n < 1 and 
the condition 



•4%V 1_/3 " =™/WPo (36) 

takes place then the intensity for the TE waves defined by formulae (|34p goes to infinity. However, 
under these conditions the absorption in the medium (and also in the walls of the waveguide) 
becomes important and the imaginary part of the dielectric permittivity should be taken into 
account. Formulae (1331) and (1341) are valid under the condition 



f" 
e o 



2 2 



,•(*) c / + c 2 

Jm, ra<- / 



(37) 



where £q = £n + ^ £ c 
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Instead of k z we can introduce the angular variable the values #m'n f° r which are related 



with the quantities km',n by the formula 



(o . a) _ muj cos i)} m . 



c 



(a,a) ' 



1 - COS ^m,' n 



(38) 



The quantities Lo m (k z ) and Aq are expressed via i?m'r? by the formulae 



n \ muJ ° fon\ 

miz) = l-frcog^ ' (39) 

Ao = ^^v^sin^- (40) 



The possible values "dm'tn are determined by formulae (fT6"|) and (|3"5j) : 



/3||6t ) „ 2 ±,/l + 6^ ) n 2 (/?J- 

'^'n ) = " -• (41) 

i + pfM 



Note that the singularity in the radiation intensity (I34p noted above corresponds to the values 
of the angular variable determined by the condition 

= ± , ± = arccos /Jy . (42) 

In the reference frame moving along the direction of the axis 2 with the velocity vu the angle 
corresponding to is equal to 7r/2. From (|4ip the following relations can be seen 

< ^ < <?o < < t, (43) 

where 

\ tf c for 0\\ > 1, v ; 

and $ c = arccos(l//3||) is the Cherenkov angle related to the drift velocity vn. Hence, for /3\\ < 1 
the waves with a = + (a = — ) are those which in the reference frame moving with velocity v» 
along the waveguide axis, propagate along the positive (negative) direction of the axis z. For 
0n > 1 the waves with a = + (a = — ) propagate inside (outside) the Cherenkov cone i? = i? c . 
From formulae ([33]) . (f34"j) it follows that the number of radiated quanta does not depend on a. 
In particular, for /3u < 1, the same number of quanta is radiated inside and outside the cone 
i? = For the case Eq = 1 we could expect this result from the problem symmetry, as in 
the reference frame moving along the direction of the axis z with the velocity vn we have a 
symmetric situation under the reflection with respect to the charge rotation plane. 

Now let us consider the radiation intensity in the limiting case of large values of the waveguide 
radius, p\ — * oo. In this limit the main contribution into the radiation intensity comes from 
large values n and we can use the asymptotic formula (see, for instance, [27] ) 



m — 1 (—1) , 
n + — s~ + ^—r- ■ (45) 



'm,n 

Replacing the summation over n by the integration and introducing as a new integration variable 
the angle for the the radiation intensity one finds 

4*> « / (46) 



9 



where 

r(°) «2, ,2 



cW c^/eosin??|l-/3||cos7?| 3 m Vi-^licos^ 

d 4m _ 9 2 ^m 2 /3|sintf J/2 / mfcsiiii? \ ^ 



dtf |1 -/?||Cos^| 3 m Vl-/?||Cosi? / 

In this limit the frequency for the radiation along given direction $ is determined by the expres- 
sion mu)o/\l — P\\ cos i? | . The expressions (|47p and ([48 p coincide with the corresponding formulae 
for the radiation in a homogeneous medium with dielectric permittivity Eq. Note that in the 
discussed limit we have the transition fim'^ — > i? and the upper (lower) sign corresponds to the 
angular region ^ i? ^ i?q ($o ^ ^ ^ vr )) where t?o is defined by formula (|44p. 

As an additional check for formulae (j33p . (|34p we can consider the special case u)$ = for 
a fixed value po- This corresponds to a charge moving with constant velocity vn on a straight 

line p = po parallel to the waveguide axis. In this limit 6m,n — ► oo and from condition (|15p it 
follows that the radiation is present only under the Cherenkov condition /3n > 1. Taking the 
limit ujq — > 0, from formulae (f33j) . (f34"|) we see that 



/(0,«)| n - 2q2 ^ J mUm,nPo/pl) r(l,a)_>n ^ 
x m,n l^o=U o n -(0) \ ' m > n v^-v 

J m +lUm,n) 

Hence, in the limit under consideration the TM waves are radiated only. The corresponding 
frequency is given by the expression vpmj n /(piy/0^ — 1). In the limit pi — > oo, by using 
asymptotic formula (|45j) for the zeros of the Bessel function, replacing the summation over n 

_/oo 

by the integration, and using the formula Q ^m( x ) = 1/2, we can see that from (|49p the 

formula for the Cherenkov radiation intensity in a homogeneous medium is obtained. Formula 
(I49p for the radiation of a charge moving parallel to the axis of the waveguide can be found, for 
example, in [29|. 

We have carried out numerical calculations for the number of the radiated quanta per one 
period of the particle orbiting, 



' max 



*£> = E E = £EEf«' (™) 

a=± n=l Q=± n=l | w m,n | 

As it has been mentioned before, the quantity Nm'n^ does not depend on a. In figure Owe have 
plotted the dependence of on the ratio pi/po for to = 24 and v±/c 0.967 corresponding 
to the energy 2 MeV assuming that sq = 3. The graphs are given for (3u = 0.9 (left panel) 
and (3\\ = 0.7 (right panel). Note that the location of the peaks in the number of radiated 
quanta for TE waves is determined by formula (|36p . In particular, for large values n by using 
the asymptotic formula (|45p we see that the distance between the neighboring peaks is given 
by formula Tryl — j3^/(mj3±) and decreases with increasing As we see the main part of the 

radiated quanta is emitted in the form of the TE waves. Similar features take place for the 
radiation on other values of the harmonic to. For the same values of the parameters, for the 
number of the radiated quanta in the homogeneous medium evaluated from formulae (|47p . 
(148|) . one finds 

2 2 

Afi « 2.683^, 9.514^, for /?„ = 0.9, (51) 

2 2 

~ 2.082^-, « 4.676^-, for 0„ = 0.7. (52) 
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Of course, in this case the result does not depend on p\. 




1.02 1.04 1.06 1.08 1.1 1.05 1.1 1.15 1.2 



Figure 3: The number of quanta emitted on the harmonic m = 24 in the form of the TE and 
TM waves per circulating period of the charge, multiplied by hc/q 2 , (hc/q 2 )Nm\ versus the 
ratio pi/po for /3|| = 0.9 (left panel) and /3u = 0.7 (right panel). The velocity of the transverse 
motion corresponds to the energy 2 MeV and the dielectric permittivity is taken eo = 3. 

In addition to the total number of radiated quanta for a given m, it is of interest to consider 
the corresponding spectral distribution. On the left panel of figureE]we have plotted the quantity 
N^'n as a function of n for pi/po = 1.02, = 0.9, for the same values of the parameters 

corresponding to figure El In this case we have n^lx = 19. On the right panel of figure H] we 
have given the corresponding frequencies for TM and TE waves. 

Now we consider an important case of relativistic charge motion in direction of the waveguide 
axis when the velocity of orthogonal motion is non-relativistic, v± <C c. This type of motion 
is realized in helical undulators. The corresponding magnetic field in Cartesian coordinates is 
given by H u = H u (— sin(fc u z), cos(fc u z), 0), where k u = 2ir/X u and A u is the undulator period 
length. The corresponding parameters for the particle orbit are related to the particle energy £ 
and to the undulator characteristics by the formulae 



v± K u v\\ 1 + K 2 K n c 

- uj = Kv\\, po = —r , (53) 



c 



7 c y 7 Z jk n vu 



where 7 = £/moc 2 , with mo being the particle mass. In formulae (|53|). K u = (q/moc 2 )H u /k u 
is the so-called undulator parameter. For example, for the helical undulator of the Stanford 
free electron laser H u = 0.23 T, A u = 3.3 cm and the electron energy £ = 43.5 MeV. For these 
values of the parameters we have K u « 0.71. As it is seen from (f33]h (j34"J) . in the presence of the 
medium the factor 1 — v jj /c 2 in the formulae for the radiation intensity in the empty waveguide 

is replaced by the factor 1 — /3 2 . This replacement leads to important influences on the radiation 
properties. These influences are essentially different in the cases eo < 1 and eq > 1. In the first 
case, even at very high energies the factor 1 — tends to finite limiting value and the radiation 
does not have the features typical for the radiation of an ultrarelativistic particle in vacuum. 
In contrast to this, when eq > 1, under the condition < 1 — /3jj <C 1 the properties of the 
radiation are similar to those for the radiation in vacuum from an ultrarelativistic particle even 
in the case when vu is not too close to c. In order to illustrate these features, in figure [5] we 
have presented the number of the radiated quanta of the TE modes with m = 1 as a function 
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Figure 4: The left panel presents the number of radiated quanta, (hc/q 2 )Nm'n\ as a function of 
n for the harmonic m = 24 and for pi/po = 1.02, /5ii = 0.9. The values of the other parameters 
are the same as those for figure [3j The black points correspond to the TE waves and the circles 
correspond to the TM waves. On the right panel the ratio wjn,«'/™o is plotted in dependence 
of n. The upper (lower) halves of the points correspond to a = + (a = — ). The values of the 
parameters are the same as for the left panel. 



of the undulator parameter K u for the undulator period A u = 3 cm and for the radius of the 
waveguide p\ = 0.5 cm. The full curve corresponds to the radiation from an electron with the 
energy £ = 25 MeV moving in the waveguide filled by air (eq = 1.00054) and the dashed curve is 
for the radiation from an electron of energy 100 MeV moving in the empty waveguide (eq = 1)- 
Note that in the first case we have 1 — f3» ~ 4.1 x 10~ 5 and in the second case 1 — /3m ~ 1.9 x 10~ 5 . 
The corresponding spectral distributions are presented in figure [6] for the value of the undulator 
parameter K u = 0.7. On the left panel we have plotted the quantity N^ 1 ^ as a function of n 
and on the right panel the corresponding frequencies are presented. The black points correspond 
to the radiation from an electron with the energy £ = 25 MeV moving in the waveguide filled 
by air and the circles correspond to the radiation from an electron of energy 100 MeV in the 
empty waveguide. 

For ultrarelativistic particles the main part of the radiated energy is in the spectral range lo > 
uj p with uo p being the plasma frequency, where the dielectric permittivity is well approximated by 
the formula £o ~ 1 — UJ p/ UJ ' 2 '• I n this regime the radiation frequencies are of the order loqj 2 and 
for the undulator parameter X u ~ 1 cm we have 7 > 10 3 . In this range the perfect conductor 
boundary condition on the waveguide walls are no longer valid and for frequencies larger than 
the corresponding plasma frequency the waveguide becomes transparent. The corresponding 
radiation intensity propagating in the exterior region is strongly directed in the forward direction 
and is described by the formulae given in [22]. As regards the frequency range u lv p the 
corresponding radiation intensities are described by formulae (|33p . (|34p . where for materials 
with £0 — 1 ^> 7~ 2 we can substitute I — (3 2 ~ 1 — Eq. By taking into account formula (j53H for 

the radius of the helical orbit, from formulae (j33H . (134)) we can see that for jm,n "C 7 the number 
of the radiated quanta Nmff is suppressed by the factor j~ 2 for both TM and TE waves, 
whereas the corresponding frequencies are practically independent on the particle energy. For 
the modes with jm]n ^ 7 we have |wm,r?| ~ cjm]n / \A — V^oPi and Nm'n^ ~ 7 -20 ". In this case 
the radiation on the TM modes dominates. 
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Figure 5: The number of the radiated quanta of the TE modes with m = 1, {hc/q 2 )N^\ as a 
function of the undulator parameter K u for the undulator period A u = 3 cm and for the radius 
of the waveguide p\ = 0.5 cm. The full curve corresponds to the radiation from an electron with 
the energy £ = 25 MeV moving in the waveguide filled by air and the dashed curve is for the 
radiation from an electron of energy 100 MeV moving in the empty waveguide. 
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Figure 6: The left panel presents the number of the radiated quanta, (hc/q 2 )N[ ^ , as a function 
of n for K n = 0.7. The values of the other parameters are the same as those for figure On 
the right panel the corresponding frequencies, lO -13 ^ 1 ^ /2tt, are plotted versus n. The black 
points correspond to the radiation from an electron with the energy £ = 25 MeV moving in the 
waveguide filled by air and the circles correspond to the radiation from an electron of energy 
100 MeV in the empty waveguide. 
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In the discussion above we have considered the radiation emitted by a single particle moving 
along a prescribed trajectory. From the point of view of practical applications the generalization 
of the obtained results for the case of the radiation from an electron beam is a next important 
step. In helical undulators the beam is bunched and the features of the radiation critically 
depend on the ratio of the bunch length to the wavelength of the emitted radiation. If the 
wavelength is smaller than the bunch length the particles in the bunch are not phase correlated 
and the total power radiated by the bunch is the sum of single particle parts. In the opposite 
limit, called the coherent spontaneous radiation regime, the wavelength of the radiation exceeds 
the length of the bunch and the particles radiate coherently. In this case the radiated power is 
enhanced by the factor of number of particles in the bunch relative to the incoherent radiation at 
the same wavelength. Here it is important to take into account that long wavelength radiation 
is suppressed by the waveguide cutoff condition and this leads to the constraints for the bunch 
length in order to have coherent radiation in the waveguide. Once the fields are evaluated, 
another interesting question is related to the influence of the radiation field on the motion of the 
radiating particles. The interaction of the beam with its own radiation can induce an additional 
microbunching with the possibility of coherent radiation from particles of the same microbunch. 
These questions require a separate consideration and we plan to address them in the future 
work. A recent discussion of the beam dynamics in undulators can be found in |30j (see also 

Another point which deserves a separate investigation is the role of the other processes of 
the particle interaction with the medium (see also the discussion in [24J). In particular, they 
include ionization energy losses, particle bremsstrahlung in media, and the multiple scattering 
(see, for instance, [32l[33]). The relative role of these processes depends on the particle energy 
and characteristics of the medium. To our best knowledge, the previous investigations in this 
direction were mainly concerned with the case of an unbounded homogeneous medium and 
the investigation of the effects induced by the presence of the waveguide requires a separate 
consideration. However, some general features can be obtained by using the corresponding 
results for a homogeneous medium. First of all, as the above-mentioned processes are absent 
in the vacuum, the radiation discussed in the present paper is the main mechanism of the 
energy losses in sufficiently rarefied medium. Next, while the ionization losses by an electron 
rise logarithmically with energy and bremsstrahlung losses arise nearly linearly, the undulator 
radiation intensity arises quadratically and, hence, dominates at sufficiently high energies. For 
a particle moving in air the ionization losses dominate those due to the bremsstrahlung for the 
energies less than 100 MeV. By using the standard formula, it can be seen that in the example 
corresponding to figures [5] and [6] the relative energy loss is ~ 1% per meter. An interesting 
possibility to escape ionization losses in the medium was indicated in [29] (see also |34j ) . In 
this paper it was argued that a narrow empty channel along the particle trajectory in the solid 
dielectric does not affect the radiation intensity if the channel radius is less than the radiation 
wavelength. From the other side, the maximum impact parameter for ionization losses is of the 
order 6 max ~ mm[(h/moc)£ / 1 , c/u p ], where / is the mean excitation energy for the atom of the 
medium, and for the channel radius larger than 6 max ionization losses are suppressed. As in other 
processes involving multiple scattering, we expect that this effect will appear in the formula for 
the radiation intensity at given frequency in the form of the additional multiplicative factor and 
there exists a critical energy of the particle below which the multiple scattering does not effect 
the radiation. We can try to estimate this factor by using the Migdal formula. For the values 
of the parameters taken in the example above this factor leads to the decrease of the radiation 
intensity by » 0.5%. 
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4 Conclusion 

We have investigated the electromagnetic field generated by a charge moving along a helical orbit 
inside a circular waveguide with dielectric filling. This type of motion is involved in magnetic 
devices called helical undulators which are inserted into a straight sector of storage rings. The 
helical undulators are used to generate circularly polarized intense electromagnetic radiation 
in a relatively narrow bandwidth. The frequency of radiation is tunable by varying the beam 
energy and the magnetic field. In this paper we have seen that the insertion of a waveguide 
into the helical undulator provides an additional mechanism for tuning the characteristics of 
the emitted radiation by choosing the parameters of the waveguide and filling medium. The 
electric and magnetic fields are presented as the sum of two parts. The first one corresponds 
to the fields of the charge in the homogeneous medium and the second one is induced by the 
presence of the waveguide. The Fourier components of the latter are given by formulae (ED , dHJ) 5 
(jlip . (|12p . In order to extract from the total fields the parts corresponding to the radiation, 
we have investigated analytic properties of the Fourier components as functions on k z . These 
components have poles on the eigenmodes of the waveguide. We have specified the ways by 
which these poles should be circled in the integral over k z . For the cases < 1 and /3ii > 1 the 
corresponding contours are plotted in figures [T] and [2j By using the residue theorem, the radiation 
fields are presented as a superposition of the eigenmodes of the waveguide corresponding to 
the TM and TE waves. The projection of the wave vector on the waveguide axis and the 
corresponding frequency are given by formulae fjl6|> . (118p . We have derived formulae (133j) and 
(|34p for the radiation intensity emitted in the form of the TM and TE waves. Limiting cases 
are considered and features of the radiation are investigated. In particular, we have seen that 
the main part of the radiated quanta is emitted in the form of the TE waves. Applications of 
general formulae to helical undulators are given. In particular, we have demonstrated that in the 
case of filled waveguide the radiation with features characteristic for ultrarelativistic particles 
in the empty waveguide is obtained for moderately relativistic particles. The radiation emitted 
on the waveguide modes propagates inside the cylinder and the waveguide serves as a natural 
collector for the radiation. This eliminates the necessity for focusing to achieve a high-power 
spectral intensity. The geometry considered here is of interest also from the point of view of 
generation and transmitting of waves in waveguides, a subject which is of considerable practical 
importance in microwave engineering and optical fiber communications. 
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